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f"*"^ i The problem of color confinement in lattice Landau gauge QCD is briefly reviewed 

, in the light of Kugo-Ojima criterion and Gribov-Zwanziger theory. The infrared 

' properties of the Landau gauge QCD is studied by the measurement of the gluon 

propagator, ghost propagator, Kugo-Ojima parameter and the running coupling 
via /3 = 6.0, 16 4 ,24 4 ,32 4 and /3 = 6.4, 32 4 ,48 4 lattice simulation. The data are 
analyzed in the principle of minimum sensitivity(PMS) in the MOM scheme. We 
' observe the running coupling ce s (q) has a maximum about 1 and decreases near 

q = 0. The Kugo-Ojima parameter, which is expected to be 1 was about 0.8. 



1. Introduction 



$— i ' In 1978 Gribov proposed that the integration range of the gauge held should 

be restricted to the region where its Faddccv-Popov determinants are pos- 
itive and that the restriction will cancel the infrared singularity of the 
gluon propagator and results in the singularity of the ghost propagator, 
which makes the linear potential between colored sources, and the color 
confinement will be realized 1 . Almost at the same time, Kugo and Ojima 
proposed color confinement criterion based on the BRST(Becchi-Rouct- 
Stora-Tyutin) symmetry without Gribov's problem taken into account 2 . 
Kugo-Ojima two-point function is defined in the lattice simulation as 



OV - ^KV) = ±5>-*(*-«'> (tr(\^D^[A„,\ b ]) xy ) (1) 



.r.y 



where u ab (0) = 5 ab u(0) and l+u(0) = Z1/Z3 = l/Z 3 . Z 3 and Z 3 are gluon- 
and ghost- wave function renormalization factor, respectively and Z\ is the 
gluon vertex renormalization factor. Kugo expressed the ghost propagator 



1 
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G(p 2 ) in the infrared asymptotic region as 3 G(q 2 ) ~ -^r-. j-t —. ». , . 

V ; V ; q 2 {l + u(0) + O(q 2 )) 

In 1994, Zwanziger analyzed various regions of the Landau gauge in 
connection with Gribov's problem, and developed the formulation of lattice 
Landau gauge 4 . From L-periodic link variables, U x ^ € SU(n) (n = 3), we 
define gauge field A^iJJ), and then there are two possible options of Ap(U), 
{/-linear (Zwanziger's) type A x ,» = (U x ^-U%^) /2\ tra ceiess part, and log U 
(ours) 7 type U x ,n = e Ax ^. Here brief review is given below in terms of 
common language to both options of the definition. In each definition, <L4 M 
under infinitesimal gauge transformation g = e e is given as 5 A 9 = D^{U)e 
with the covariant derivative defined as D^{U)4> = S {11^)3^ + [A^,^], 
where operations <9 M and ~^ on the scalar <p to give vectors are defined as 
(dn4>)ic,n = H x + M) - <P( X ) and ( < / >A ')x,/i = (<P(x + n) + 4>{x)) /2, respec- 
tively, and the operation S{U^) on a vector is defined for each op- 
tion as S(U,)B, = (1/2) { (U X)M + U^/2, B x .,} ^ in [/-linear 

definition, S(tf M )B M = S(A X ,„)B X ^ = {{A x ^2) /th{A x ^/2)}B x ^ in logCZ 
definition with A Xjfl = adj Ax ^ = [A Xtfi , ■ ]. 

The Landau gauge 3A = can be characterized 6 such that 
5Fjj(g) = for any Sg, in use of the optimizing functions Fjj(g); 
F u(g) = E x ,„^ { 2 - + U 9 J)} in [/-linear definition, and F v (g) = 

Ex,^ tr ( A9 xJ A l,^) = ( A l\ A l) in lo S u definition. The variation 
AFjj(g), under infinitesimal gauge transformation g~ 1 6g = e, reads as 
&Fu(g) = -2(3A 9 \e) + (e| - 3D(A 9 )\e) H . The fundamental mod- 
ular region is specified by the global minimum along the gauge or- 
bits, i.e., A L = {U\ A = A(U), Fu(l) = Min g Fu(g)}, A L c Q L , 
where £Il is called as Gribov region (local minima), and FIl = {U\ — 
3D{U) > , 3A = 0}. So far, all field variables are supposed to be L- 
periodic. Zwanziger further defined the core region El as a set of the 
global minimum points of Fjj{g) in the extended gauge transformation 
g = e UJ e 6x where U and u> are L-periodic, and constant belong to 
an arbitrary Cartan subalgebra, i.e., [9^,6 V ] = 0: El = {U\ Fu(l) = 
Min g=elu e exFu(g)} C Al- Putting g = e UJ e 6x with arbitrary w and 9^ 
of order e, one has the variation of Fjj(g) up to 0(e 2 ) for U 6 Si as 
AFu(g) = 2(A fi \e fi ) + (cj'\-3D\cj') - {DO^-dD)- 1 ^) + (e^\S(U^\9^ > 
0, where u>' = uj — (—3D)~ 1 D6. The non-negative sum of the third 
and fourth terms is expressed by — (O^H^du) with the horizon ten- 
sor = —D^{—3D)~ 1 D V — 6^S(U^). Taking the trace of the operator 
H^ u with respect to the normalized constant colored vectors rf^ a = d^X a 
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with trA at A b = 5 ab , one defines the horizon function H(U) as H(U) = 

E„M' a \H w \vZ' a ) = E^jvr\- D ^- dD y lD ^K a ) - (" 2 - i ) e ( u ) = 

h(U)V where (n 2 - l)E(U) = J2 X ,^ ^{X^ S(U x ^)X a ). Thus one has for 
U G 3 L that A^ = V~ l Y, x A x,v, = and H(U) < 0, where V = L 4 . 
Zwanziger hypothesized that the dynamics on 3^ tends to that on in the 
infinite volume limit, and derived the horizon condition, statistical aver- 
age ^/i(t/)) = 0, in the infinite volume limit. Taking the Fourier transform 
of the tensor propagator of the color point source, (—D JJ/ (—dD)~ 1 D I/ ') xa b , 

one has G^( P )S ab = {e/d){ Pfi p u /p 2 )S ab - {8^ - (p^/p 2 )}u ab (p 2 ) , where 
e = (E(U))/V and dimension d = 4. Putting Kugo-Ojima parameter as 
u(0) = — c and comparing limp „->+o with (h(U)} = one finds that 
the horizon condition reduces to 

<^^>=(5) + ^- 1 > c - e =^- 1 ) ( c -5)^ (d - = °- (2) 

Kugo-Ojima's and Zwanziger's arguments emerge to be consistent with each 
other provided the lattice covariant derivative meets with the continuum 
one e/d = 1. Accordingly, Zwanziger derived independently the same char- 
acteristic singular behavior of the ghost propagator 4 as Kugo's 3 , both per- 
turbativcly in the sense that the diagrammatic expansion was used. Recent 
Dyson-Schwingcr analyses 9 and numerical simulations are to be noticed in 
this respect as well. It should be pointed out here that our numerical simu- 
lation does not pursue the selective core region dynamics, but rather checks 
the Zwanziger hypothesis in our standard method 8 . 

2. The gluon propagator 

The gluon propagator and its dressing function are defined numerically as 

d,M) = ^^J2 e ~ iq{x ~ yh <M^Mv)) = (<W - ^f)D A {q 2 ) (3) 

and Z(q 2 ) = q 2 DA(q 2 ), respectively. The infrared behavior of Z{q 2 ) is 
parametrized as (q~ 2 ) aD , and we observed that ao < for L > 24. 

As a kinematical structure of the Gribov region, Zwanziger gave a 
bound 5 on the constant mode = V~ 1 ^2 X A X ^ for A <G Ql, as 

2 7T 

\AJ < — {tan( — )} = a(L) with a lattice size L,V = L 4 , and a = -v/3/2 in 
a L 

SU(3) and a = V2 in SU{2), respectively. This fact can be derived irre- 
spective of the gauge field options. The Zwanziger's bound was checked 
to hold in case (3 = 6.4, L = 48, as \A^\ < 0.005 vs. a(L) = 0.10 
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for each sample. The propagator at zero momentum has an extra fac- 
tor V, as D A (0) ~ VA 2 . Although W(A^) = 0.23 ± 4.2, the propa- 
gator has the value 3D A (Q) = ^V(E M .„V ) ~ 150 ± 30 ( averaged 
over samples). It implies that the gluon propagator is infrared finite. 
Since there is a strong cancellation between samples, the sample average is 
fVE„W = 0.16±31. 

Here we review Zwanzigcr's strong argument for the vanishing of zero 
momentum connected Green functions in the infinite volume limit 5 . Let 
a partition function be Z(J) = J dAp(A)e JA , and for constant J x = h, 
it reads that Z{h) = J dAp(A)e h - A ^ = J dAp{A)e vh7i , where A^(0) = 
a d J2 X A^(x). Let W(h) = log Z(h), and then the connected Green function 

VG { c n) (Q,...,0) = ((2(0))")| c is given as G< ,l) (0, 0) = ^^W(h) . 

Note from positivity of p(A) that eT vh<7 ^ < Z(h) < e vh,j{L \ Putting 
w{h) = W(h)/V, we obtain |w(/i)| < ha(L) — * as L — > oo. Hence it 
proves that Gc (0, ...,0) — > in the infinite volume limit. However, it 
should be noted that this infinite volume limit is only proved directly for 
connected Green functions. 

In Figl, we plot the gluon dressing function of 24 4 , 32 4 and 48 4 lattices. 
The infrared properties of the gluon dressing function is studied in the prin- 
ciple of minimum sensitivity(PMS) 13 ' 14 in MOM scheme. The parameter 
y = 0.02227 is fixed at q = l.QIGeV and q 2 D A {q 2 ) = Z(q 2 , y) | y = .02227 is 
compared with the lattice data. 




2 4 6 8 10 12 14 



Figure 1 . The gluon dressing function as 
the function of the momentum q(GeV). 
/3 = 6.0, 24 4 ,32 4 and f3 = 6.4, 48 4 in 
log U version. The fitted line is that of 
PMS in MOhl method. 




0.5 1 1.5 2 2.5 3 



Figure 2. The ghost propagator as the 
function of the momentum q(GeV). (3 = 
6.0, 24 4 ,32 4 and (3 = 6.4, 48 4 in logf/ 
version. The fitted line is that of PMS in 
MOM method. 
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3. The ghost propagator 

In use of color source |A a x) normalized as tr(A a x|A h a;o) = 5 ab 5 XyXo , the 
ghost propagator is given by the Fourier transform of 

D G b (x,y) = (tr^xli-dDiU)}- 1 ^)) (4) 

where the outmost () denotes an average over samples U. The ghost dress- 
ing function is defined as G ab (q 2 ) = q 2 DQ ab (q 2 ), and in the infrared region, 
it is parametrized as (p~ 2 ) aG . Although the ghost dressing function is 
monotonic in the case of L < 32, the lowest momentum point of L = 48 is 
suppressed. 

We compare in Fig. 2 the result of the PMS method in MOM scheme 14 
D G (q 2 ) = -Z g (q 2 ,y)/q 2 \ y=0 .o2i42 and the lattice data of /? = 6.0, L — 24,32 
and P = 6.4, L = 48. 

4. The Kugo-Ojima parameter 

We directly measured the Kugo-Ojima parameter on 16 4 ,24 4 ,32 4 and 48 4 
lattices. Results arc summarized in Tablcl. 



Table 1. Lattice size dependence of the Kugo-Ojima parameter c, the trace e 
divided by the dimension d and the deviation parameter from the horizon condition 
h. The suffix 1 corresponds to the (7-linear definition, and 2 to the log U definition. 






L 


Cl 


ei/d 


hi 


C2 




h 2 


6.0 


16 


0.576(79) 


0.860(1) 


-0.28 


0.628(94) 


0.943(1) 


-0.32 


6.0 


24 


0.695(63) 


0.861(1) 


-0.17 


0.774(76) 


0.944(1) 


-0.17 


6.0 


32 


0.706(39) 


0.862(1) 


-0.15 


0.777(46) 


0.944(1) 


-0.16 


6.4 


32 


0.650(39) 


0.883(1) 


-0.23 


0.700(42) 


0.953(1) 


-0.25 


6.4 


48 








0.720(49) 


0.982(1) 


-0.26 



In Fig3, the value c is plotted as a function of log Z3(1.97GeV) for 
j3 = 6.0, L = 16,24 and [/-linear and logU definitions of the gauge field, 
and for j3 — 6.4, L = 32,48 and logU definition. The crossing point of 
the extrapolation of the points corresponding to the log U and U— linear 
definitions of the gauge field suggests c in the continuum limit c = \ — Z\ /Z3 
is consistent to 1. 

5. The QCD running coupling 

Using the renormalized gauge field A^(x) = Z3 -1 ' 2 A"(x) and the ghost 
field c b (x) = Z~ 1/2 c b (x), one can show that the vertex rcnormalization fac- 
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Figure 3. The Kugo-Ojima parameter 
c as the function of log Zz(l.97GeV). 
/3 = 6.0, 16 4 , 24 4 in [/-linear (diamond) in 
log [/(triangle) and j3 = 6.4, 32 4 ,48 4 in 
log U version (star). 



Figure 4. The running coupling ot a (q) of 
= 6.0, 24 4 (box), 32 4 (triangle), /3 = 6.4, 
32 4 (diamond) and 48 4 (star) as a function 
of momentum q(GeV) and the result of 
the PMS method in MOM scheme. 



tor of the ghost field Z\ = Z g Z 3 Z3 = 1, implies g 2 Z^Z 3 is renormalization 
group invariant. 

We estimated the running coupling from the two-point function (gluon- 
and ghost-dressing functions) 

« S (P 2 ) = {g 2 /(4n)}Z(p 2 )G(p 2 ) 2 ~ (pa)- 2(aD+2ao) - (5) 
Dyson-Schwinger approach predicts 9 ac = 0.5953, an = —1.1906 and 
a s (0) = 2.972, while the MOM scheme approach of Orsay group and the 
PMS method 14 predict that a s (0) is small and close to 0. 

Table 2. The exponent of gluon dressing function near 
the zero momentum old, near the pa = 1 a' D , the expo- 
nent of ghost dressing function near the zero momentum 
olq and od + 2a<3 in log U type simulation. 





L 


OLD 




a.G 


an + 2«g 


6.0 


32 


-0.375 


0.302 


0.174 


-0.03(10) 


6.4 


48 


-0.273 


0.288 


0.109 


-0.06(10) 



ALPHA collaboration 12 derived scheme independent running coupling 
in the Schrodinger functional method. Our data of a s for q > lGeV is 
consistent with those of Schrodinger functional. 



6. Discussion and Conclusion 

In the lattice simulation of MOM scheme 10 , the running coupling decreases 
in the infrared region and Ajj-g — 295±20MeV was reported. Orsay group 
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claims that the non-perturbative gluon condensates reduces the A-p^- in the 

M OM scheme from 295MeV to 238MeV and becomes consistent with that 
in the Schrodinger functional method 11 . 

In the PMS 13 one assumes that there is an infrared fixed point and one 
bridges the gap between ultraviolet asymptotic free and the infrared non- 
perturbative regions. In our analysis of lattice data we considered optimum 
coupling constant y at q = 1.97GeV that corresponds to 1/a of j3 = 6.0 and 
observed that the PMS in MOM scheme fits the lattice data, q > 0.5GeV 
in the ghost propagator, and q > lGeV in the gluon propagator and the 
running coupling. 

In the infrared region, we also analyzed the data via the contour- 
improved perturbation series which is expressed by the Lambert W 
function 15 . We find that the result is close to that of the static quark 
potential derived by using the perturbative gluon condensate dynamics ay 
shown in the analysis of hypothetical r lepton decay coupling constant 16 . 
Results will be presented elsewhere. 

This work is supported by the KEK supercomputing project No. 02-82. 
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